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THE INWARD BULGE TYPE BUCKLING OF MONOCOQUE CYLINDEES 
Y - REVISED STRAIN ENERGY THEORY WHICH ASSUMES 
A MORE GENERAL DEFLECTED SHAPE AT BUCKLING 


By N. J. Hoff, Bertram Klein, 
and Bruno A. Boley 


SUMMARY 


A strain energy theory is developed for the calculation of the .crit- 
ical load for the inward "bulge type of general instability of reinforced 
monocoque cylinders subjected to pure bending. The deflected shape at 
buc kling is assumed to be represented by an expression containing eight 
free parameters in addition to the two characterizing the. wave lengths in 
the circumferential and axial directions. The theory is applied to two 
representative cylinders of the GALCIT test series and to two of the 
PLBAL series. The critical stresses calculated are 8.3 to 22*9 percent 
higher than the experimental values. 


INTRODUCTION 


When a reinforced aluminum-alloy monocoque fuselage is subjected to 
bending moments, such as those caused by the aerodyna m ic loads acting 
upon the tail surfaces, structural failure is likely to occur by the 
buckling of the stringers on the compression side of the fuselage. The 
half wave length of the deflected shape is eq u al to the spacing of the 
ring fiames when the frames are comparatively rigid and are spaced far 
apart. The critical stress can be increased by spacing the rings closer, 
but at the same time the cross-sectional dimensions of the rings must be 
made smaller in order to keep the weight of the structure unchanged. 

Such a trend can be noticed in the development of the modern monocoque 
fuselage from 1930 until the present time. When, however, the distance 
between rings and the cross section of the ring decreases beyond a certain 
limit, the rings no longer have sufficient bending rigidity to act as 
rigid supports for the stringers. In such a case failure occurs by general 
ins tability in which several stringers and rings are involved simultaneously. 

The general instability of a reinforced thin-walled curved shell was 
first investigated by Uschou in 1935 (reference l) -vfho, at the ggesticn 
of H. Wagner, worked out and applied to monocoques the buckling leory cf 
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orthotropic shells established in 1932 by Flugge (reference 2) . A 
similar theory -was worked out independently by Taylor (reference 3 ) 
in 1935 * In. these investigations a uniformly distributed compressive 
lo ading was assumed. It was believed that the results would also apply 
in the case of bending provided that the wave length in the circumfer- 
ential direction, were smal 3 enough so that the variation in the intensify 
of the compressive stress caused by the bending moment could be disre- 
garded. At the same time the wave length must be greater than the spacing 
of the reinforcing elements. Otherwise the structure cannot be considered 
as an orthotropic shell. The results of these theories were compared 
later with experiments carried out at GALCIT in 1939 (references 4 to 8). 

It was found that the buckling loads predicted were several times as higi 
as those observed. 

Another form of funeral instability, caused by bending, was investi- 
gated by Heck in 1937 (reference 9) whose work is a generalization and 
application to reinforced monocoqu.es of an earlier theoretical analysis 
of JBrazier in 1927 (reference 10 ). This type of instability consists of 
a flattening of the rings and is of no practical importance * in aiuminum- 
alloy reinforced monocoque fuselages . ' 

The problem of general instability due to bending -was approached from 
another angLe by Hoff in 193 8 (reference H) • The most highly compressed 
stringers were considered as columns elastically supported by the rings, 
and the buckling load was calculated by the strain energy method. The 
deflected shape was assumed to be mainly an inward bulge symmetric with 
respect to the most highly compressed stringer. The stringers were 
assumed to be distributed uniformly around the circumference but the rings 
were treated individually. The buckling load was minimized with respect 
to two parameters which were the lengths of the wave in the axial and the 
circumferential directions. Tests carried out with two cylinders Bhowed 
that the inward bulge type of general instability was possible and that 
it occurred approximately at the stress predicted by theory. A non- 
dimensional structural parameter A was found, the magnitude of which 
determined whether failure would occur by general or panel instability. 

The extended experimental investigations undertaken at GALCIT begin- 
ning in 1938 proved Hoff's theory to be too conservative for practical use. 
In GALCIT 'a theoretical development Yon Harman and Tsien (reference 12) 
found that the classical linear theory had to be replaced by a no n l ine ar 
theory in order to predict correctly the buckling loads of nonreinf orced 
cylinders. Because of the complexity of the nonlinear theory even in the 
case of nonreinf orced cylinders the idea of finding a theoretical solution 
of the buckling problem of reinforced cylinders was given up at GALCIT. 

An empirical formula based on dimensional analysis was established which 
was in good agreement with the results of the GALCIT tests. 

In 19^3 Hoff .(reference 13 ) published a new version of his theory. A 
review of his earlier work disclosed that there the circumferential wave 
length, not expressed explicitly, was always equal to the total circum- 
ference of the cylinder after the buckling load was minimized. On the 
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other hand, in the tests the inward "bulge was restricted to the most 
highly compressed regions of the cylinder. In the new publication the 
discrepancy "between experiment and theory was attributed to the fact that 
in the original work the shear strain energy stored in the sheet covering 
was neglected. On the other hand, the strain energy of shear was found to 
"be unduly large when it was assumed that the expressions describing the 
distortions of the stringers and rings were valid also for the sheet. 
Apparently there are patterns of distortion for the sheet which involve 
less strain energy, particularly when the sheet is already in a buckled 
state at the moment when general instability occurs. 

The difficulties involved in finding the actual distortions of the 
.sheet covering were circumvented through a semiempirical approach. The 
wave-length parameter n was assumed so as to obtain an agreement between 
the buckling load observed in the GALCIT tests and -the predictions of 
the theory. In reality the wave length must depend upon the Variation cf 
the shearing rigidity of the panels of sheet around ihe circumference of 
the c ylin der, and the shearing rigidity is certainly influenced by the 
compressive and shearing stresses prevailing in the panels. It was con- 
cluded, therefore, that n was dependant upon two parameters charac- 
terizing the state of the sheet, namely r/d » nfl s r/t. Tha connection 

between these two parameters and n was established from an evaluation of 
the GAIX3IT test results, and the correctness of the dsbumptions was sub- 
stantiated by the fact that a family of smooth curves was obtained. "With 
the aid of these curves the buckling load could be predicted in a simple 
manner. 

The purely theoretical analysis was resumed in 19^5 by Hoff and main 
(reference lh) at PIBAL. It was found that the shear strain energy stored 
in the sheet covering was not unreasonably large when it was calculated on 
the basis of an average shear angle defined by the relative displacements 
of the four corners of each panel. Such a treatment disregards details of 
the states of stress and, strain in the panel w/nfl should be based on experi- 
mental values of the shearing rigidity of panels subjected to arbitrary 
compressive loads. In th$ a&aence of experimental data a connection between 
shearing rigidity and compressive load was assumed which was hoped to 
represent fairly the actual conditions. Other improvements upon the 
original treatment of the subject were the replacement of the expression 
def inin g the shape of the Vave in the 1 circumferential direction by a new 
expression which satisfied more rigorous boundary conditions the 
replacement of the integration around the circumference by a summation. 

The new results were in good agreement with the GALCIT test data. 

An experimental investigation of the shearing rigidity of reinforced 
mono co que cylinders carried out at PIBA1 by Hoff and Boley (reference 15) 
in 19^6 , however, resulted in the measurement of -considerably higher 
shearing rigidities than expected in the case of curved panels subjected to 
compressive loads exceeding the critical loads. The use of such high v alue s 
of the shearing rigidity in the revised theory just discussed would h&ve 
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resulted In "buckling loads considerably larger th a n those observed. At 
this stage of the development the authors began to entertain same doubts 
whether the problem could be solved by a linear theory. 

Their conviction that a linear theory would predict correctly the ( 
buckling load in general instability had been based on the following 
observations. The inward bulge proper was always restricted to a region 
of the cylinder where the sheet covering was in a buckled state when 
general instability occurred. Because of the increased movability of the 
wavy curved sheet the re inf orced monocoque cylinder was reduced essen- 
tially to a grid consisting of stringers and, rings connected somewhat 
loosely by the buckled sheet. The question arises now whether the support 
rendered by the rings to the stringers has a linear characteristic. In 
experiments carried out at PIBAL with-' reinforced' monocoque cylinders a. 
linear connection was found between a load applied radially to a ring and 
the deflection of the point of attack. Under such a load the entire ring 
distorted, whereas in tests undertaken at GALCIT in order to clarify some 
aspects of the n online ar theory the nonlinear support was furnished by a 
semicircular element, the ends of which were rigidly fixed. In the 
general instability tests both at GALCIT and PIBAL the rings were observed 
to distort • slightly around the entire circumference so that the end points 
of the main bulge were definitely not fixed. Because of these .consider- 
ations one more attempt was made at PIBAL to explain the phenomenon of the 
inward bulge type of general instability by a linear theory. The results 
of this investigation are presented in the present report. 

The strain energy method is known to yield too high buckling loads 
if the assumptions regarding the deflected shape deviate, frpm the actual 
shape of deflections. An effort was made therefore to make the new 
assumptions more elastic than were the earlier ones by the use of a greater 
number of parameters that can adjust themselves in the process of minimi- 
zation. Altogether nine free parameters were incorporated in the expression 
describing the deflected shape, one of which is indeterminate as in all 
buckling problems. To the .remaining eight parameters may be added the tiro 
parameters defining the wave length in the circumferential and axial 
directions . 

The strain energy and the work of the external forces were calculated 
by summation as in the earlier publication (reference 1^) , and in the 
deter min ate on of the shear strain energy stored in the sbeet covering the 
experimentally obtained values of the shearing rigidity were used. The 
large number of parameters did not permit the development of an expres- 
sion in closed form for the buckling load but four representative cylinders 
were investigated numerically and -the calculations resulted in buckling 
stresses in' satisfactory agreement with the test results. The authors 
believe that the agreement obtained between theory and experiment 
substantiates the claims for an essentially linear character of the prob- 
lem in the form it was set up in these calculations . Details of the 
behavior of the buckled panels, however, cannot be explained in all 
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probability by a linear theory. It is quite likely that better agreement 
between theory and experiment can be obtained if an even more complicated, 
or else a better, shape is assumed for the deflections. 

The calculations presented here are too lengthy for use in structural 
design. The writers believe that the .simple formulas developed in 
reference 13 are those most suitable for practical calculations. The 
value of the parameter n should be taken from a diagram published by 
Hoff, Boley, and Nardo (reference 1 6 ), which is a revised version of the 
diagram presented in reference 13 . The differences between the two dia- 
grams are the incorporation of new test data obtained in a series of tests 
at PIBAIi and the 'correction of the earlier test results for the effects 
of stresses exceeding the elastic limit. The importance of the theory 
presented in this report is twofold. First, it can be used in the calcu- 
lation of n values in that region of the parameters r/d and e TnR - ir r/t 

which is not easily accessible to experimental methods. In such a mariner 
a more comprehensive n diagram can be constructed than the one available 
at present without recourse to a very extensive test program involving 
cylinders of the size of the fuselages of modem transport airp lane s. 

Second, it establishes the fundamental soundness of the linear approach 
upon which the method suggested for practical calculations and the n diagram 
are based. 

For the sake of completeness other publications dealing with general 
instability are now listed. In 1937 Nissen (reference 17) compared test 
results obtained with reinforced cylindrical panels and complete cylinders 
subjected to compression with.Dschou's theory. Satisfactory agreement -was 
found after an empirical factor was incorporated. Byder (reference 18 ) 
attempted in 193 8 to improve Taylor's theory by the use of eupirical factces. 
In 19^2 Tsun Kuei Wang (reference 19 ) developed a theory of the flattening, 
of monocoque cylinders subjected to compression, but the applicability cf 
his results is Impaired by the fact that he neglected the shear strain 
energy stored in the sheet covering. At PIBAL additional theoretical 
investigations related to the inward bulge type Instability were carried 
out by Hoff and Klein (reference 20) in 1944 and a series of cylinders ms 
tested by Hoff, Fuchs, and Cirillo (reference 21) in the same year. A 
new general instability theory was worked out by Wang (reference 22) 
in 1946 in which the wave form was described by infinite series. The shear 
strain energy stored in the sheet was not calculated but its effect was 
taken into account by means of an empirical multiplying factor depending 
upon the size of the panels. 

For his contribution to the calculations presented in the present 
report the authors are indebted to Merven W* Mandel. This investigation 
was conducted at the Polytechnic Institute of Brooklyn under the sponsor- 
ship and with the financial assistance of the National Advisory Committee 
for Aeronautics. 
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SYMBOLS 


a Q Indeterminate coefficient 

- a/b,c,d,e,f ,g,h parametric coefficients 


A 

A str 

d 

E 

®red 

®sec 

®tan 

Gr 


cross-sectional area of stringer plus effective width, of 
sheet 

o 

area of stringer 

width of panel measured along circumference 

Young's modulus 

reduced modulus 

secant modulus 

tangent modulus 

shear modulus 


Gq shear modulus of sheet covering at zero compressive load 

\ 

G eff effective shear modulus 

h centroidal height of stringer 

i index denoting position along circumference 

I moment of inertia 


I r moment of inertia of ring cross section arid its effective 

width of sheet for handing in its own plane 


I str ro 


moment of inertia of stringer cross section for tending in 
radial direction (about its tangential principal axis) 


Ia-bj. moment of inertia of stringer cross section and effective 

width of curved sheet for tending in radial direction 
(about its tangential principal axis) 

Ig-fcr moment of inertia of stringer cross section for tending in 

"k° tang ential direction (about its radial principal axis) 


Istr-fc moment of inertia of stringer cross section and effective 

width of curved sheet for tending in tangential direction 
(about its radial principal axis) 
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J 

L 

m + 1 


M 


n 

Pl> P2 
P 


Q 

r 

s 

S 

t 

U 




U ( 


ah 


U 


str. 


J atr 4 


2w 

w n 

w r 

v t 

W 

x 


index denoting position along axial direction 
length of wave in axial direction 
distance "between adjacent rings 

0 

number of ring fields involved in failure 

function appearing in strain energy of bending of rings 

parameter defining wave length in circumferential direction 

polynomial functions of a,b,c,d,e,f ,g, and h 

force acting an stringer and its effective curved sheet at 
buckling 

function appearing in shear strain energy 
radius of cylinder 

number of stringer fields involved in bulge 
total number of stringers in cylinder 
thickness of sheet covering 
strain energy 

bend in g strain energy stored in rings 

shear strain energy stored in sheet 

radial bending strain energy stored in stringers 

I 

tangential bending strain energy stored in stringers 
effective width of curved sheet 
rotation of tangent of ring 

radial displacement of a point on a ring or a stringer 
tangential displacement of a point on a ring or a stringer 

O 

work done by applied forces at buckling 
axial coordinate 

coefficients used in calculation of shear strain in a 

panel due to rotations and displacements of its corners 


“r* “t 
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P 

7 

8 


cr 

€ ° r sh 

e cu 

*ft 

A 

<P 

<D 




angle subtended "by a panel 
shear strain 

shift of neutral axis from horizontal diameter at buckling, 
percent of radius 

normal strain in a stringer plus its effective curved sheet 
at buckling 

•maximum compressive strain at buckling 

compressive buckling strain of sheet panel 

buckling strain of a nonreinf orced circular cylinder under 
uniform axial compression 

buckling strain of a flat panel of sheet in compression 
structural buckling index 
angular coordinate 

function of cp appearing in expression for w r 
function of <p appearing in expression for 

DE7SL0PMEHT OF A REVISED BUCKfflTG THEORY WHICH INCLUDES 
EIGHT ARBITRARY PARAMETERS 


Revised Deflection Pattern 


The shape of the bulge at buckling is determined mainly by the 
radial deflections. The following expression is chosen to represent the 
radial deflections w r occurring at buckling: 

w r * a o[® r ^ a ^) sin 2 («x/L) + ® r (c,d,e) sin 6 (nx/L) 

‘ + ® r (f ,g,h) sin 10 (jtx/L)j (l) 

provided that 0 

0 <cp < (ic/n) 

(la) 


0 < x <, L 
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where 


$ r (c,d,e) = [c cos op + d cos 2 ncp + e cos 3 111 ? 

+ ( 2 . 666 . ..c - 2 * 333 ...d + 1 . 777 ***e) cos 4nq> 

+ ( 1 . 666 . . .c — 1 * 333 * • • d. + 0 . 777 . . . e) cos 5 ^^ (lb) 

The value of $ r (l,a/b). is obtained by replacing c, d, and e in equation (ib) 
byl, a, and b, respectively; $ r (f,g,h) is obtained by replacing c, d, and e 
byf, g, and. h, respectively. 

Also v r = 0 when cp > it / n and/or x >L (lc) 

The notation and the sign conventions are shown in figure 1. 

The deformations of the rings are assumed to be inext ens i onal . The 
condition for inextens ionality is: 


w r « - dv-t/dcp 


( 2 ) 


Equations (l) and (2) determine the tangential deflections w^: 

Wj. = - (a^n) J®j.(l, a,b) sin 2 (rtx/L) + $ t (c,d,e) sin^(itx/L) 

+ 4> t (f,g,h) sin 10 (rtx/L)[ ( 3 ) 


provided that 


0 < cp < (tf/n) 

. ' - T'OC* Z.S. 

0 < X <L 

.‘r xl.- r..i L Ul'> -- '• 


. .o? .-.x 7 

.6 s. reli’.''? - id" c 


' : (3a) 
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■where 

<&.j.(c,d,e) = jjj sin ncp + (l/ 2 ) 4 sin 2 rap + (1/3) e sin 3119 

+ (l/ 4 ) (2.666* . 1 c — 2.333. . • ! 1.777. ..e) sin 4 imp 
+ (1/5) ( 1 . 666 ... c - I.333 ... 4 + O.777... e) sin 5 n<p] ( 3 b) 

The value of* $.j.(l,a,b) is obtained by replacing c, 4 , an 4 e in equation (Jb) 

by 1 , a, an 4 b, respectivelyj $.j.(f,g,h) is obtaine 4 by replacing c, 4 , an 4 e 
by f , g, and h, respectively. 

Also Wj. = 0 when cp > it/n an 4 /or x > L (3c) 

A few explanatory remarks are now given regar ding the choice of the 
4 eflecte 4 shape. In orAer that the bulge be symmetric about the bottom 
stringer cp = 0 the following conditions must be satisfied there: 

1 . The tangential displacements must vanish: 

w^ = 0 when «P = 0 , for all values of x ( 4 a) 

and 

2 . The radial deflection pattern must exhibit a horizontal tangent: 

(dw r /kq>) = 0 when 9 = 0, for all values of x (4b) 

Further, in order that there be a smooth transition between the bulge and 
the nondistorted part of the cylinder at 9 = (it/n) , the following conditions 
must be satisfied there: 

3 * The tangential displacements must vanish: 


( 4 c) 


w t “ 0 


when 9 = (rt/n) , for values of 


x 
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4 . The radial displacements must vanish: 


w r = 0 when 9 = (it/n) , for all TalQ.es of r (4d) 


5 * There must "be no sudden change in the direction of the tangent: 



when 9 


(it/n), for al l values of x ( 4 e) 


and 


6. There must “be no sudden change in the curvature: 


(d%/&p 2 ) = 0 when 9 = (rt/n), for all values of x (4f ) 


Finally, in order that there he a smooth transition between the bulge 
and the nondistorted part of the cylinder at x = 0 and x = L, the fol- 
lowing conditions must be satisfied there: 

7 « The tangential displacements must vanish: 


w^ a 0 when x = 0 (or L), for all values of 9 (4g) 

8. The radial displacements must vanish: 


L), for all values of 9 (4^) 


in the direction of the tangent: 

(or L) , for all values of 9 ( 4 i) 


and (3) satisfy all the boundary 
to ( 4 i) . (The arbitrary function 
of x which would normally appear in equation (3) as a result of the inte- 
gration of equation (2) is zero as a consequence of co ndi tion (4a).) 
Typical examples of the deflection patterns of a stringer and a .r ing , taken 
from the results of the calculations, are shown in figures 2(a) and 2(b). 


w r = 0 


when 


x = 0 (or 


and 


9 * There must be no sudden c hang e 


(dv r /d:^ = 0 when x = 0 


It can be shown that equations (l) 
conditions enumerated in equations (4a) 
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Calculation of Strain Baergy 


Str»^ energy stored, in rings .- The strain energy stored in any one 
ring Is 

nit/n 

U = (1/2) |(El) r /r3j I [w r + 2 &p (5) 

ll -Jt/n 

4 

If the value of w r is substituted from, equation (l) and the strain energy 
is s umme d over all the rings, the following expression is obtained: 



(l,a/b) + 9 ^- 



(El) r sin^^yi^ 


' J=1 


(EC) r sin 2 (rtxy^ 

J=1 


+ |$ r (f,g,h) + ^’(f^h)] ^ 


(H) r sin^^xy^ ^ 


2 


dcp 


(6) 


where m is the total number of rings included in the wave length L. 

The integration yields a result in closed form. The summations contained 
in equation (6) can he evaluated algebraically as functions of (m + l) . 
They are given in table I. IT all the rings have the same bending 
rigidity (El) r , the total strain energy U r stored in all the rings becomes 


TT r . (V2)(#*3) (S0 r * MtW> X + M( °' 4 ’ e) ^ 






( 7 ) 



I 


viie re 


M(o,d,e) = jfito 2 - l) 2 + (2.333 ...) 2 (l6n 2 - l) 2 + (I.333. . . ) 2 (fen 2 - j) 2 ] d 2 

- 2 [(2.333...) (I.777...) (l6n 2 - ]) 2 + ( 1 - 333 -..) < 0 . 777 ...) (fen 2 - l) 2 J de 
+ [(fe 2 -2) 2 + ( 1 . 777 ...) 2 (l6n 2 -j) 2 + (0 .777 . . .) 2 (25n 2 - l) 2 ]e 2 

- 2 1 ( 2 . 666 . . . ) (2.333. . . ) ^.6n 2 - l) 2 + ( 1 - 666 ...) ( 1 - 333 -.. ) (fen 2 - a) 2 ] cd 

+ 2 j(2. 666...) (1.777...) (l6n 2 - 3) 2 + (1.666.. . ) (0.777.. •) (fen 2 - ]) j ce 

+ [(n 2 - l) 2 + (2.666...) 2 (fen 2 - 2) 2 + (1.666. .. ) 2 (fen 2 - }) 2 J o 2 (7a) 


The value of M(l,a,b) is obtained by replacing c, d, and e by 1 , a, and b, respectively* M(f ,g,h) Is 
■ obtained by replacing c, d, and e by f, g, and h, respectively* is obtained by replacing d 2 

by ad, de by (1/2) (ae + bd), e 2 by be, cd by (l/2)(d + ac), ce by (1/2) (e + bo), and o 2 by 05 

M [J' d O ia °' btained V replacing d 2 by dg, de by (l/2)(dh + eg), e 2 by eh, od by (l/2)(cg + df), 

/ a a TVS 

ce by (l/ 2 ) (ch + ef), and c 2 by cf. ^Similar considerations can be used to find nQ* * 


k 


,h 

00 


1 
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Strain energy stored In stringers .- The strain energy stored, in the stringers "because of "bend in g 
in the radial direction is 

s 


<W r - Y_ (1/2) (m) etr r 
1 





cbc 



(B) 


where the s umma tion is extended over all the stringers contained in 'the "bulge. Substitution and 
integration yield 



0|)» ri (l,a,li)» ri (f,g,li) + (^^)» ri (o,d, 0 )* ri (f,g,H) 


(9) 


The strain energy stored in the stringers "because of bending In the tangential direction Is 


<W t -X 

1 



( 10 ) 
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Substitution and integration yield 


<W t .- (1/2) ^C a W 2 X ( ” )3tr t 

* J " 


(^§§ji) + (¥)V 1 ' a ’ 1,) V°' 4 ' e) 

/ 

(^D *t t Cr,8Pi) + (i§g) « tl (<=.!-») ® tl (f,e,i>)] 


(11) 


Because both the * functions a-nd (KE) 0tr and (BI) s -t r . fc vary from stringer to stringer, the summations 
appearing in equations (9) and (11) have to be evaluated numerically. 

Strain energy of shear stored in sheet . - The shear strain energy per unit volume is taken as 
the average effective shear modulus G 0ff ’multiplied by one -half the square of the average shear 

strain y in the panel. The latter is calculated from the displacements of the four corners of the 
panel. Then the total strain energy of shear stored in the sheet is 

- 1 sh - d/2>H7 2 <W V 1 (12) 


where Irjtd is the volume of a panel. 

The effective shear modulus depends upon the geometric and mechanical properties of and the 
average normal strain in the panel. Its value was taken from the empirical curves established 
earlier at PIBAI and presented in figure 24 of reference 10. 
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The average angle of shear 7 was calculated. from the equation 

r = (i/l,) [°*(- + v ri+i;J + v rijJ+1 - w r±+1 ^ 

+ Kl (- w %,j - + + w t 1+1)J+ i) 

+ | (- w nij3 - »n 1+1> j + "iii ; j +1 + 'Vi + i, j + i) j to) 


where the first subscript refers to the circumferential location and. the 
second, to the axial location of the corner of the panel, as shown in 
figure 3. The rotation w n of the tangent of the ring is given by the 
relation 


w n = (l/r) (Sw r /dq>| (14) 

The values of the factors 04., 04., and. were calculated, from the 
equations 


a r = p(b.l + 0.000238p 2 ) 

04. = - 0.5 - o. 01666... p 2 ) 
a n = ~ (l/P) (o«5 + a -(J 


( 15 ) 



Substitutions yield 


m .|_i 


U sh = ( X /2) ^(td/l^Go ) 2^- ^ eff /G^ j[sin 2 ^L_ . gin 2 * ^^[ |%(l,a,b) 


[f^rih: - sl “ 6 2 ^ J rr i ] t! i (o ' 4 ' e) + ^“srfr " stal ° ~ ^ t r ) ] a i (f,8,1>) l' 


(l£) 


"Where ^ is the number of stringer fields involved ip. one-half of the symmetric bulge. The meaning 

nf* ■KH© amnrfhnl Ch f r* A a\ 4 a 
— — \- w J'*.? w / **- r 


Qi(c,d,e) 


°«r£- W 0 ' 4 '®) + ®ri+l (c ' d ' e) ] + (K|/ n ) 



+ |<%[(d/r)n - (1/n) t ri '(c,d,e) - (l/n)® ri+i '(c,d,e)~ 


(16a) 


It is possible but cumbersome to sum up the .trigonometric functions describing the deflected shape 
in the axial direction contained in equation (l6). It was found more convenient to carry out these 
summations numerically. The same practice was followed as regards the 4> functions. 
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Work Done by External Forces 

Equal a-nfl opposite forces are assumed, to "be acting at the x = 0 
a-nri x = L ends of each stringer. The distance "between the points of 
application of these forces shortens at buckling. The product of 
this short ening and the force is the work done by the force. The total 
exte rnal work is the sum of the work done by the forces acting on each 
stringer contained in the bulge: 


W = 




(&w r /ax) 2 + (8*^4 2 J 


dx 


( 17 ) 


Tha a^r! ai strain distribution is assumed to be linear. Then the 

external force acting upon the ith stringer at buckling, may be expressed 
in terms of the buckling strain e cr of the most highly compressed stringer 

and is: 


*i 



6 

or 


( 17 a) 


where rB is the shift of the neutral axis at buckling. Substitutions and 
integration yield: 


I 



\ 



1 

1 


+ (ll) [®ri (1 ' a ' l) ®r 1 (r ’ e ’ h) + 



“*1 


(18) 


j The summations encountered in equation (18) were performed numerically. 

| Calculation of Buckling Load 

The Buckling condition is 

u r + u str r + ^str^ + u ah * w 

t 

i 


(15) 


£ 


os. £Z YSN& 
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•where the Tallies of the strain energy quantities and of the work must "be 
taken from, equations (7), (9), (ll) } (l6), and (l8). Equation (19) was 
solved for sqj., which is a multiplying factor in the expression for ¥, 

and the result minimized by means of the following procedure. 

Integral values of s and (m + 1) first were chosen for the 
number of stringer and ring fields included in the bulge. On the basis 
of ttese tentative values the M, and Q functions reduced to quadratic 
expressions of the eight parameters a to h. Next € cr was assumed) this 

permitted the calculation of the. shift of the neutral axis, the effective 
widths of sheet (taken to act with the stringers), and the moments of 
inertia of the stringers and made it possible to read the values 
of (Geff/Go) *cam the appropriate graph. All necessary summations were 

then carried out. Substitution of all the results in equation ( 19 ) made 
it possible to obtain a solution for e cr in the form: 


S = ^1 (ljS-jb,c,d,e,f ,g,h) 

Cr £2 (l,a,b,c,d,e,f,g,h) 


( 20 ) 


.lore p x p 2 u. qpatortio poljnandals to the sigHt para^ters a 

to h. Minimizing this expression for with respect to each of these 

eight parameters is equivalent to setting 


_ Pi _ dpi/fta _ dpi/8h _ _ dpi/dg _ dpi/dh 

€ cr pg ^pg/Sa dp^/db dpg/dg dp 2 /dh 


whore the partial differential coefficients of p^_ and p g are linear 

functions of the eight parameters. Equation (21) represents nine connections 
between e cr and the eight parameters. They were solved by a trial-and- 

error procedure. First a value of 6 was assumed, and a ket of values 

cr 

of the eight parameters was dete rmin ed by solving eight linear equations. 
These values were substituted into the original quadratic expression to 
obtain a calculated value for 6 Qr . This procedure was repeated with new 

assumptions for until the calculated, value was reasonably close to 

the one assumed. In order to locate the absolute minimum value of e cr t it 

is necessary to perform the calculations for a number of different choices 
of s and (m + l) • An .example of the numerical calculations is shown in 
the appendix. 
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COMPARISON OF THEOET AND EXPERIMENT 


The theory was applied to four monocoque cylinders tested earlier 
at GALCIT (references 4 to 8 ) and PTBA1 (references 1 6 and 21). 

Dimensions of the cylinders showing rings and stringers are given in 
figure 4. Typical "buckling patterns obtained in the calculations are 
shown in figure 2 , and a tabulation of the numerical results is presented 
in table II. Theoretical and experimental strains rather than bending 
moments are compared because of the following reasoning. It was found in 
the calculations that the buckling strain remains practically constant if 
a normal strain distribution differing materially from linearity rather 
than a linear normal strain distribution is assumed to- exist at buckling 
This is so since the- most highly compressed stringer (cp = 0) contributes 
mainly to the bending and work terms and also since the other strain 
energies are little affected by changes in the strain distribution. 

Although the critical strain is little affected by changes in strain distri- 
bution, the value of the bending moment may vary considerably. In order 
to calculate the latter, it is thus necessary to know the strain profile 
existing at buckling. Such knowledge was not available for the GALCIT 
cylinders. It was found convenient to assume a linear strain distributim 
'for the PIBAL cylinders as well a p for these. 

The axial wave length predicted by theory for PIBAL cylinder 11 was 
larger than the total length of the test specimen. It was found- experi- 
.mentally in references 4 to 8 that if the length of a cylinder is at least 
twice its diameter its buckling load is independent of its length. This 
is not quite true for PIBAL cylinder 11 but the effect of the smaller length 
should be unimportant . 

The four cylinders were so chosen as to cover the widest possible 
ra n ge in physical properties and to obtain critical stresses below or not 
far in excess of the proportional limit. The latter requirement was set 
up- in order to prevent the variation in the effective modulus from 
attaining major importance. As may be seen from table II the theoretical 
buckling strains are greater than the corresponding experimental values as 
they should be when the strain energy method is used. The deviation is 
small; it amounts to 22.9 percent of the experimental value in the worst 
case and to 8*3 percent in the best case. 


CONCLUSION'S 


The strain energy theory developed for the Calculation of the critical- 
load for the inward bulge type of general instability was found to give 
critical stresses 8.3 to 22-9 percent higher than the experimental values 
obtained with four representative cylinders tested at GALCIT and PIBAL. 
This agreement is considered satisfactory for most practical applications 
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Taut it is likely that "better agreement could "be obtained by wa wta of more 
comprehensive, or more suitable, assumptions for the deflected shape at 
buckling. The numerical calculations involved are too lengthy for use in 
structural design. 

lor this reason it is suggested that practical, calculations be carried 
out with the aid of the procedure suggested in the paper entitled "General 
Instability of Monocoque Cylinders” by H. J. Hoff (Jour. Aero. Sci., vol. 3D, 
no. 4 , April 1943, PP* 105-114, I 30 ) and in NA.C/LTH Ho. 1499- The sig- 
nificance of the present investigation is twofold, lirst it permits the 
calculation of values of the parameter . n, needed, in the practical pro- 
cedure suggested, in a region where experimental data are not available and 
too expensive to obtain. Moreover the satisfactory agreement between 
results of the present theory and experiments Indicates • that the general 
instability phenomenon is essentially a linear problem can be calculated 
by means of a linear theory. 


Polytechnic Institute of Brooklyn 
Brooklyn, H. Y., March 3 , 1947 
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APPENDIX 

CAIX3ULAIEEQNS FOE GALCIT TEST CYLINDER ,6? 


In this appendix details of the calculations performed in determining 
the critical strain for GALCIT cylinder 67 are shown. The geometric and 
mechanical properties for this cylinder are: 


1 ^ y in. . • 
Tj in. . . , 
S .... a 1 
t y in. . . 

Astrj 

I rj ^ •, 

Ip,-*-,. , in. 

s-CTro^ ^ 

^trto^ 111 ‘ 
Eqj psi • • 

G 0, P si • • 


...... 4 

10 

12 

. . . . 0.012 
. . . 0.0324 

219.3 9 x 10-7 
. 374 x 10“ 6 
. 563 X 10-6 
. 10.5 X 10° 
. . 3.9 X 10° 


' I 

The reported experimental critical strain is 29*5 X 10 . At the 

outset of the calculations e cr is taken as 33 x 10“^. The value of 

16 is calculated to he 1.692 inches on the basis of a linear strain 
distribution. This. permits the setting up of table A. 



TABES A 


£ 


[jer ^ 33 X 10"^j rB « 1.69s dLnf] 


(1) 

(2) 

( 3 ) 

(«0 

( 5 ) 

(6) 

( 7 ) 

(S) 

. ( 9 ) 

(10) 

(n) 

(12) 

1 

< 

2 w 

x atr r 

x otr-t 



^eff 

VfC^cr) 

e boo 

®tan 

3 J?ed 

0 

1 

s 

33.OO X 20-4 

2 Q.22 

3I.89 ■ 

1-257 

i. 3 te 

1.721 

9.996 x 10-^ 

6 . 08 S 

6.366 

—— 

15 

13.28 

8-55 

O.3975 

■3575 

•305 

0.04497 

.0458s 

.04961 

0.04497 

.04057 

.oeQ 4 

10 

10.9 

10.9 

6.69 

9.9O 

IO.5 

8.25 

9.98 

IO.5 

RJPjQFQJj 


TABLE B 

&-g 


1 

oos(2rfi/4) 

oos(4iti/4) 

cos(6rti/4) 

00s (&ti/4 j 

OOB (lOffi/4) 

0 

1 

1 

1 

1 

1 

1 

0 

-1 

0 

1 

0 

J 

Multipliers 

c 

4 

e 

2.666.. .0 

1.666. • >o 

■f/w A f A /I a\ 
— 




J.555. . .^a.1 .'rrr. ..a 

’"•—'III — 

4/) m YT / a a > A 
— * — 

1 

sin (2*1/4) 

sln(4«i/4) 

sin(6*i/4) 

Bin(8jti/4) 

sin(10«i/4) 

0 

0 

0 

0 

0 

0 

1 • 

1 

0 

-1 

0 

1 

Multipliers 

0 

0.5<1 

0*333* * 

0-666...C 

0.333. **o 

for *,.(0,4,0) 




-0.5833. . .440.444. • .e 

■0 .266* • *4+o *155* • *® 

Multipliers 

0 

24 

3* 

20*666* • «o 

8.333- --o 

tor -O r '(o,l,e)/n 




-9*333* • *4+7 • ill* • • 0 

-6.666. . .44-3*888. . .e 
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Column (l) refers to the stringer station. Column ( 2 ) is the strain at 
the locations of the stringers. Column ( 3 ) is the effective width of 
curved sheet calculated from the equation: 


2v = (1/e )(d/r)|o.3t + 1.535 [(t/d)(er - 0 . 3 ^ 1 ^ (Al) 


Columns (4) and ( 5 ) give the moments of inertia of the stringers plus 
their effective width of curved sheet calculated from the equations: 


L str, 


" I str ro + 


fr * (t/2) : P*° a H} 2 , 0 . a [(gif ) 2 /gitj 2 gift m 


(l/2wt) + (1/A) 


•^str^ ~ ^str-tQ + (l/l2). (2w)^t 


(A3) 


Column ( 6 ) indicates the ratio of the actual strain in a panel when the 
monocoque "buckles to the "buckling strain of a panel of sheet in compres- 
sion calculated from Eedshav 1 s formula: 


«cr oh * (W 2 ) + /(W 2 ) 2 + s cu 2 = 2-2 x 1°" 4 W 


The values of column ( 7 ) are read from figure 2k of reference 15 . The 
entries in column ( 8 ) are given "by Agtr + 2wt. Those in column ( 9 ) can 
"be computed from columns (2) and ( 8 ). If the strain in some members is 
above the proportional limit it is necessary to use reduced moduli of 
elasticity. In calculating the bending of the stringers the Yon TfrCrm^n 
reduced modulus was used; in calculating the work of the external 
forces', the secant modulus was used. These values are based on the 
curves of reference 23 and given in columns (12) and (10) of table A. 

For the evaluation of the 4 functions it is convenient to set up 
a tabular arrangement . For s = if-, it is given by table B. By UBing 
this table the polynominal ® r (c,d,e) for 1=0 is found by multiplying 

in each column the expressions in the row below the first double line 
by the numbers in the corresponding columns listed in the first row of 
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the table and adding the results of all the products* In a s imilar -marmoT* 
the polynomial for ® r (c,d,e) for i = 1 is obtained as -well as the two 

expressions for each $t(c,d,e) and -® r * (c,d,e)/n. The results are pre- 
sented in table C. 

The function [j-$ r ^(c,d,e) + ®ri + i(c,d,e) and the two others needed 
to calculate ^(c^d^e) can be determined with the aid of table C by simply 
subtracting or adding the polynomials in adjacent rows for each c olumn of 
that table. In doing this it must be remembered that all the ® functions 
’va n ish when 1 = 2 - The results are not shown hare . "When S = 12, s = 4, 
and n = (12/4) = 3 equation ( 15 ) gives: 


0 * = 0 . 05239 ^ 

|o^|/n = O .165905 

. |a^j (d/r)n = O.OO 6854 


(A5) 


The value of (c,d,e) determined from equation (l 6 a) is recorded in 
table D. The other ® and Q functions are derived by replacing c,d,e in 
tables A to D by l,a,b, and f,g,h, respectively. In the theory the # and Q 
functions appear squared. It is thus necessary to square the polynomials 
of tables C and D. The results are not included here. The M functions 
are taken from table IH j in this instance the ones corresponding to n = 3 
are needed. 

The next step in the calculations is to assume a value of (m + l), 
the number of ring fields involved in failure. The value of (m + l) was 
taken as 7 and the summations appearing in equations ( 6 ) nnfl ( 16 ) were 
evaluated numerically. (Table I afforded a check for the former.) These 
results as well as the coefficients of the $ functions for the r stringer 
bending (equations ( 9 ) and ( 11 )) and the external work {equation (l 8 )) a re 
tabulated in table 17 . These numbers are then multiplied by the appropriate 
one of the M, $2, and Q , 2 functions to yield for each strain energy and 
for each value of i a quadratic form in the eight parameters consisting 
of 45 terms. These are to be found in table Y. For the shear strain energy 
a constant value of G e ff /^0 O ’ 3575 is ta ken so that the quadratic expres- 
sions for i = 0 and i = 1 can be combined into one quadratic expression 
given in table Y. 

For purposes of calculation it is convenient to multiply each of equa- 
tions ( 7 ), (9), (11), (16), and (l 8 ) by ' — , where E se ~ is 

(l/ 2 )ao 2 (* 2 /L)E 8eCo 
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the secant modulus of the most highly compressed stringer, and solve for 
]> 

6 cr X 10 instead of e cr . In order to calculate pj_ and p 2 in 

equation (20), each row in table Y must be multiplied by its appropriate 
factor. These are given as follows with numerical values calculated from 
data previously presented*. 












TABLE I 

Matrix 1 far x lcA - 35 


a 

b 

0 

4 

e 

f , 

e 

b 

-EBB 

692.2146826 

-464.4565336 

588.7416155 

497.5350394 

-332.9098263 

-4&}. #773353 

409.2661605 

-273.4394178 

-818.5987754 


343.7051960 

385.9414027 

-332.9099263 

246.3400096 

316.480^30 

-273.4394178 

•202.2678 766 

540,99094®. 



566.3080636 

-475.9672414 

305.0229474 

491.3603230 

-412.7507501 

262.4374082 

700.4413908 




404.3549033 

-267.3909489 

-412.7507501 

351*6459540 

-231.4902753 

-588.7416157 





198.3617120 

262.4374082 

-231.4902753 

172.0563584 

385*9414027 






444.5869972 

-372.6801402 

235.4-15&705 

575.5740447 







319.5626629 

-209.158&705 

-4^.8773352 








I56.I653638 

316.4808330 









974.5640600 

Values of the 

unfrnoyna 








0.620032608 

-O.899355851 

-0.018836025 

2.015640161 

3.011252608 

0.571555092 

--0.457014134 

-1.547918125 



■4to lawir part of tbs matrix baa Loan emitted because of symmetry. 
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The expression € cr = is obtained by multiplying idle fore- 

going factors by the entries in table Y term-by-term and adding. A 
first guess for e x 10^ equal to 35 is made. The corresponding 

eight equations are contained in table E. They were solved by Cr out * s 
method (reference 24). The resulting values of the eight parameters 

were substituted into equation (20) to yield e cr = = 36.79 X 10”^ . 

On using 34-5 as a second guess for x 10^,. a value of '34.67 
was obtained which was considered close enough to the new guess to termi- 
nate this procedure. 


A graphic method may be used to determine accurately the critical 
strain after several steps in the approximate solution of the buckling 
equations have been undertaken. At buckling the value of -the determinant 
formed by the coefficients of the eight equations with the right-hand-side 
elements transferred to the left-hand side must vanish* If a plot of the 
value of the determinants against the strains assumed at each stage is 
made, then the critical strain will correspond to the zero value of the 
determinant of this curve (extrapolated if necessary) . Also the value of 
the determinant will be positive or negative, the sign depending on 
whether the assumed strain is below or above the critical strain. Justi- 
fication for these statements may.be found in reference 25* In this case 
such a procedure yielded ^ x lCr = 34.4 (see table II). 


Instead of* repeating the calculations for new values of (m + l), it 
was found more suitable to adopt the following scheme in minimizing the 
buckling strain in the axial direction. It can be shown that if 
<t(l,a,b) /<t(c,d,e) and. $(l,a,b)/$(f ,g,h) are independent of cp, then the 
ring, stringer, and shear energies divided by the external work beoome 

nearly proportional to (m + l) 2 , l/(m + l) 2 , and (m + l) 2 (l - cos Y 

' m + 1/ 


respectively* Since this requirement is approximately fulfilled, it is 
believed that the critical strain found by proportioning the three strain 
energies according to the foregoing factors will be close to the critical 
strain derived in the rigorous manner. This procedure is permissible 
provided that the critical value of (m + l) calculated in this manner is 
near the value for which the calculations were carried out rigorously, 

7 in this instance. This method entails the assumption that the number 
of ring fields involved in failure is not an integral number. The trans- 
formation from (m + l) = 7 to (m + l) =7.63, the approximate critical 
value obtained, is indicated as follows!. 
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Strain 

Energy 

(m + 1) = 7 

(m + 1 ) = 7.63 

Ring 

207.75^ 

247.020 

Stringer, radial 
1 = 0 

216.256 

182.018 

i = 1 

67.735 

57.011 

Stringer, tangential 
i = 1 

10.003 

8.40$ 

Shear 

i = 0 and i = 1 

216.226 

218.734 

1-1 

717.974 , 

713.202 k 

e cr x 10 ^ = 3^*67 

€ cr x 10 * = 34.44 

P 2 

20.7Q9 

20.709 


Since the value e c ~ X 1(A = 34.4 Instead of 34.67 corresponds to crit- 
ical strain for (m + 1 ) = 7 , the correct critical strain for (m + l) = 7.63 
may he found by the following equation: 

€ cr x 1C,k m ( 34 . 44 / 34 . 67 ) (34.4) = 34.17 (A 6 ) 


As seen from table II, assuming s = 2 or s = 6 ( (m + l) =7 in. each case) 
gives values of the critical strain higher than the foregoing one obtained. 
Consequently 34.17 x 10”^ may be taken as the critical stra in . 
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TABLE II 


XABOLAXICff QF BSSQUB 


Cylinder 


(m + 1) 

Strain used 
at cutset 

Strain used 
In solving 
equations 

Sign at 
deter- 

Trrtnfln-h 

Magnitude ocf 

flfltmdpaat 

Strain, 

2s 

P2 

Value 

from 

graph, 

« 

ox* 

Critical 
(»+ 1) 

Heir 

*«p 

Eaperl- 

aaotal 

(a + 1) 

Experi- 

mental 

'OX - 


\ '««® / 

PIBAL 


■ 


fit x 10" * 

D 

E.309 X 10 12 

RgHK 

- 






cylinder 

6 

fl 

a 5 x 10-* 

25 

■ 

l.Oflt x 10 32 








U 


■ 


25.6 

■ 

0.438 X 10 10 


26.1 X 10'^ 

9.55 

22.59 x lo-* 

6 

20.86 x 10-* 

8.3 

PXBA 1 




19 

- 

4.923 X 10 32 

20.98 





• 


oylindsr 


fl 

17.9 

10.7 

- 

1 . 5 t 3 x up 2 

18.79 

18.6 

7.95 

18,13 

7 

16 

13-3 

54 


B 


18 

+ 

8:812 x 10 7 

22.48 




- 







lfi.O 

B 

12.943 X 10 12 

18.04 

17.2 

7.26 

17.2 

8 

14 

22.9 

cylinder 

it 

7 

19 

16. V 

9 

8.316 X 30 32 

IB.32 







30 


| 


16.7 

■ 

5.445 x jo 32 

18.74 






V 

QALCIT 


■ 


35 

+ 


98.6 




B 



cylinder 

9 

fl 

35 

35 

- 

2.771 X 10 32 

38-79 




■ 



67 




34.5 

- 

5.749 x 10 32 

3^-67 

34.4 

7.63 

3*. 17 


29.5 

15. 8 





to 

- 

-170.65 x io 12 

60.86 




9 
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TABUS HI - M X 10" 1 * 


n 

ad,ag,dg 

at),d©,gh 
1/2(09 + dh) 
l/2(ae + M) 
l/2(ah + Tjg) 

* V ,* 2 

"be/bh,eh 

a # cd,fg 
l/2(cg + df) 
1/2 (ao + d) 
l/2(af + S) 

1 b,ce,fh 
l/2(oli + ef ) 
l/2("bc + e) 
1/2 Cbf + h) 

i 2 * 2 - p 2 
1 ,0 

o,f ,cf 

2 

1.9629 

-2.6628 

0.9849 

-4.6476 

3.1521 

2.7729' 

2-5 

3.8714 

-5.2521 

1.944365 

-9.16362 

6.216645 

5.467005 

2.666. . . 

4.7083643 

-6.38755378 

2.36517478 

-H.14384979 

7.56o48l32 

6.6483359 

3 

6.725340741 

-9-123990123 

3.379410699 

-15.916029629 

10.799091357 

9.495229630 

3*333*** 

9.24650181 

-12.54446474 

4.64728517 

-21.88093464 

14.84723739 

13.05364357 

4 

I6.O2545 

-21.74153333 

8.05667222 

-37.919000 

25.73191666 

22.62125 

5 

31.3757200 

-42.56746666 

15.7775644 

-74.23455999 

50.37909333 

44.285440 

6 

54.28867036 

-73.65386173 

27.30298313 

-128.44081481 

87.16917901 

76.62231482 
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TABLE XV 

COEOTICIMTS OF M. *2. AND «2 RMCTI0N8 FOR 
STRAIN ENERGY AND EXTERNAL WORK 

[(» + 1) ■ 7 j 


Ring 

strain energy 


Stringer 
strain energy 


Shear 

strain energy 


External 

work 


Bln 4 (rtj/7) 


sln 12 (rtj/ 7 ) 


H Bln^Jtj/T) 

2 Z ain 8 (rtj/7) 

2 |I sin“(*j/r) 
2 .21 sin 1 ” (*J/ 7 ) 


Coefficient of G^fl^a,!)) 
Coefficient of 4 2 (c,d,e) 
Coefficient of * 2 (f,g,h) 
Coefficient of *(1^,1) *(o,d,e) 
Coefficient of GCl^a,!)) <f(f,g,h) 
Coefficient of ®(o,d,e) *(f,g,h) 


5Z Bln 2 (jtj/7) - sln2jt(j + lj/xjf 

sln^(jtj/7) - sln®it(,J + 1 )/tJ £ 

_sln 12 (nj/7) - sin 12 *(j + l)/7] 2 
2 X jjih 2 (irj/7) - sin^J + 1 )/tQ eln 6 (*j/7) - sln 6 *(j + l)/7j 

2 21 {ain 2 (itj/7) - sln^fj + l)/7] ein^fitj/^) - sin 10 ir(J + l)/f 

2 H ^in 6 (xj/7) - sin 6 *^ + l)/T] [sln 10 (nj/ 7 ) - Bin 10 *^ + 1)/X 


Coefficient of 
Coefficient of 
Coefficient of 
Coefficient of 
Coefficient of' 
Coefficient of 


t|(l,a,t) + (l/n)3t 2 (l/a,t) 

4>§(o,d,e) + (l/n)^ t (c,d,e) 
v^(f,s,h) + (l/’n)% t (r ,s,h) 

® 7 (l,a,h) «> 7 (c,d,e) + (l/n) 2 # t (l,a/b) ® t (c J1 d # e) 
• 7 (l,a,l>) » (f,g,h) + (l/n) 2 * t (l,fc,b) * t (f,g,h) 


$ 7 (o,d,e) ® 7 {f,g,h) + (l/n) ® t (o,d,e) 


2.625 

1.579101562 

1.218159722 

3.8281249952 

3.158203122 

2.7424926744 


2 

6.690625 

I4.I8304443 

3*75 

3.28125 

17.88574218 


O.658892846 

0.992954708 

O.967843396 
1 . 2354s4o88 
1.080996077 
1.816155626 


0-5 

0.73828125 

0 . 92735^05 

0.9375 

0.8203125 

1.571044921 


a 


jSACA 













00 


lAEtB-T 

onADBmc topm nr a,b,c,<i,e,r b 


(s - 12* e - 4} (» + 1) - 7] 


Strain energy 

1 

Oonatant 

a 2 

b 2 

o 2 

d 2 

e 2 

f 2 


Bing 

Btringer, radial | 

Stringer, tangential 
. Shear 

0 

1 

1 

24.924977779 

56.888888860 

14.222222222 
3 • W555536 
•.0278825156 

17.654010445 

14.002022202 

22.222202202 

0.142202222 

.0150011066 

8.870953065 

25.2839506 

6.32098766 

.063209B8 

.0178996484 

14.99393194 

106 

40 

12.25 

.<1377873637 

IO.619996069 

- 

76.5623 

0.49 

.020330026I 

3.336432733 

67.111111111 

21.777777778 

0.217777778 

.0242362317 

11.366696791 
403 .4288193 
100 .8572048 
23 .214301212 
.040936443 


Wort: 

0 

1 

14.222222222 

3.6543e099 

3.5555?5556 

5.559506175 

6.320987695 

1,582)02745 

21 

5.395833337 

li ^8958337 

9*333333333 

2.335925928 

26.37803818 

6.777690372 




e 2 

h 2 

ah 

2ao 

ad 

, 2ae 

2af 


Bing 

Stringer, radial £ 

Stringer, tangential 
Shear 

0 

1 

1 

8.102532482 

100.8572048 

15T.3Q338B3 

1.008572048 

.020035027!) 

4.116658618 
L7?. 3016975 
44.8254044 
0.448254268 
.0260926766 

.23.950474073 

-37.925925925 

-23-703703703 

0.189629629 

.0061666862 

-60.928930848 

-IO6.666666667 

-66.666666667 

-2.666666667 

-01443IQI32 

25.745444991 

26.666666667 

41.666666667 
0.266666667 

.02812707.49 

-34.927774644 

-71.111111111 

.44.444444444 

0.335555538, 

.0115625366 

-50.266054467 

-93-333333333 

-58.333333333 

-2.333333333 

-.0126271366 


Work 

0 

1 

' 6.5945095 
10.31124840 

11.723572556^ 

2.9341496826 

-9.48l4GQ.46l 

-5.920658435 

-26.666666667 

-16.740740738 

6.666666667 

10.424074078 

-17.777777778 

-11.101234566 

:S:S« 




ag 

Hah 

2bo 

2bd 

be 

£bf 

£bg 

bb 

Bing 

Stringer, radial j 

Stringer, ' tangential 
Shear 

I 

21.239S92125 

23-333333333 

36.45^33333 

0.233333333 

.•004*111905 

-28.815414092 

41.340271559 

iko.PwePiWP- 

-34,927774644 
-71 .nnnnnn 

12.936806563 

47.407407407 

H.85185I85I 
O.II8318525 
• 0333618409 

34.105724038 

124.444444444 

3I.IUIIIIII 

-1.555553356 

-.0733039792 

-28.EQ.5414092 

-62.222222222 

-38.888888889 

O.jllllimill 

10.67286542 

41.481401461 

10-37037037 

0.103703703 

.0293666108 


35.55555556 

-1.777777778 

-0837759762 

44.444444444 

O.355535536 

.QH5625366 

.0101172106 

.0113172196 

Work 

0 

1 

ism 


35.555535556 

8.839506169 

-17.7Y7777778 

-U.I01234566 


3l.nnnill 

7.734567898 


10-37037037 

2.595473254 




SoGt ’Off m vam 











I 


qiHIWin BOBt a *,>,0,4,11,7,6, JHD t - OmoWU 
[a ■ a - *j + 


Strain mutter 

“ 

od 

c* 

Of 

*8 

Boh 


ttf 


Btria^r, r*il*l }► 

BtertaoBO 1 , t«n|]iiitTii 1 
SbMT 

9 

- 85 . 15302791 ^ 

- 196.0 

- 122.3 

4.9 

-.010430662 


26.040597698 

Sr’Sg 

A 

43-649594657 
-306.73 _ 

.-.02121*530 

^g.6l£42i&9 

BO 

g§j 


¥Qtk 

B 

-Q 1 '0 

- 13.103333333 

es.o 

S.osmiiTi 

4*16875 
U. 462204866 

-^■^38*83 

VMM 

- 1 *^ 

- 8 . 7*2222222 

— 

- 





Eerf 

Sot 

oh 

fa 

■C^H 

Bh 

Uag 

Strlaanr, radial ^ 
BMagsr, traawtlta 

fflD0 

a 

1 

1 

1B.4441ST7M 

ll&TS&S 

Uli 

Btefl 

BB1 

9 . 26800906 * 

pofi.nmiin 

36.357778 

.flipytMJift 

|g( 

B 

-U.ll** 6 Bl *8 

- 268 . 9323*6186 

■^®7 

•OCom&m 

Vox* 

B 

wf** 

^ 5.751666666 

^ 13 . 60^07127 


-S».79iW6tf6 

-nlB.eo3fii»OTfi7 

lo.mmim 

4*97^79*735 

EE3 


-IO .963079598 



• 

1 

0 

a 

So 

f 

% 

Eh 

Sins 

trtriago, radial ”1 

Mrine«r, tngraMM, 
Bbrar 

B 

Sp!3 

BPS] 

0*000000000 

•036679717 

m 

Jl*5tatTl^7 

itt.aa2iaa 

m 

m 

iS:**a5*5 

5i.miniii 

-l’33B3i»335 

-073303979 

--007696540 

-M&oyzi 


lforfc 

B 

-lV.tS 22 S 22 fia 

4.92639306 

18 . 9669*966 

e£ flaaoaetff 
KS. OOOOOOOOO 

£.8318131892 

- 26 . 66666666 $ 

- 16 . 7407*074 

’HESS 

23.333333333 
3 ’993370370 


«.nmiin 

V.TUtO^B 

7.734367693 




U) 

vo 
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SECTION A-A 

STRINGERS NOT SHOWN 
OUTSIDE OF BULGE 


FIGURE I.- THEORETICAL DEFLECTED SHAPE 


HACA 03J No. 1505 



X 


X-L 


X-yL 


(Q) BOTTOM STRINGER. <P = 0. 



(b) RING AT X=fL. 


FIGURE 2r TYPICAL DEFLECTIQN PATTERNS. GALCIT 
CYLINDER 67. S- 12, s«4,(m + |) = 7. 



FIGURE 3.- DEFORMATION OF CORNERS OF A PANEL SHOWING 

NOTATION USED IN CALCULATING SHEAR STRAIN. 
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GALCIT CYLINDER 67 


FIGURE 4r MONOCOQUE CYLINDERS. 





